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1. Introduction 



The world-sheet description of the non-critical string theories is given by the two- 
dimensional quantum gravity (for reviews see flHl]). Considered as a theory of coupled 
Liouville and matter fields, the 2D QG has much larger symmetry than each of its two 
components. Turning on the gravitational field erases all the information about distances 
and instead of correlation functions one obtains 'correlation numbers'. One of the manifes- 
tations of this larger symmetry of 2D QG is the so called 'ground ring' integrable structure 
discovered in the early 90's It has been shown by Witten that for c = 1 matter 

coupled to Liouville theory, there exists a ring of ghost-number zero operators, with respect 
to which the local observables, represented by ghost-number one vertex operators, form a 
module. Kutasov, Martinec and Seiberg ^ have pointed out that the action of ground 
ring leads to to a set of recurrence equations for the correlation functions of the primary 
fields (the closed string 'tachyons'). The ground ring structure is particularly interesting 
also by the fact ||,^ that it resembles the integrable structures observed in the matrix 
models of 2D QgI. 

Bershadsky and Kutasov [|ll| extended the ground ring structure to the case of 2D 
QG with boundary, or c = 1 open string theory. In addition to the two bulk ground ring 
generators, it involves two boundary operators, which generate the 'boundary ground ring'. 



In 1 11], the action of the boundary ground ring is used to find recurrence relations for the 
correlation functions of boundary operators, or open string 'tachyons'. The compatibility 
of these (over determined) relations allows to evaluate all structure constants and all cor- 
relation functions, and to reproduce elegantly the results obtained previously by Coulomb 
gas integrals | |T^ - |TB| . 

Originally the ring relations were applied only for the so called 'resonant' amplitudes, 
which are the only non-vanishing ones in a theory with vanishing bulk and boundary cos- 
mological constants. In the full theory, the resonant amplitudes give the residues of the 
'on-mass-sheir poles of the exact amplitudes. It was conjectured that the perturbations 
like the cosmological term are described by certain deformations of the ground ring struc- 
ture, but the exact field-theoretical meaning of this conjecture remained obscure. Only 



after the impressive developments in Liouville CFT from the last several years [17-21] it 
has been realized that the ground ring structure is an exact symmetry of the full CFT 
constants [^. Indeed, the conformal bootstrap used to derive the exact expressions for 
the bulk and boundary structure constants in p!9| - |21[] is essentially the same procedure as 
the one used to establish the ground ring relations. The only technical difference is that 
while the identities obtained in pure Liouville theory follow from the truncated OPE with 



^ Very recently, new evidences about the deep connection between these two integrable struc- 



tures were put forward in the interesting paper by Seiberg and Shih [10]. 
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degenerate Liouville fields, in 2D QG one considers products of degenerate Liouville and 
matter fields. 

In this paper we consider the boundary ground ring structure in Liouville quantum 
gravity with non-vanishing boundary cosmological constant, or equivalently, in a 2D string 
theory in presence of a FZZT brane. Extending the the technique developed by Bershad- 
sky and Kutasov we obtain a set of functional recurrence equations for the n-point 
boundary correlation functions. Crucial role in our analysis will play an observation made 
by V. Fateev, A. Zamolodchikov and Al. Zamolodchikov |19|| about the nature of the 
degenerate boundary Liouville fields, which we sketch belowB. 

In boundary Liouville CFT there are two cosmological constants // and hb, the bulk 
and the boundary one. ^b/^^- The observables are meromorphic functions of /x^, with 
a branch point singularity at /i^ = — /i^, where fi^ ~ ^/JI is minus the boundary entropy 
produced by the fluctuations of the Liouville field in the bulk. The branch point can be 
resolved by introducing a uniformization parameter s 

= A^B cosh(7r6s) (1-1) 



where b is the Liouville coupling constant. FZZ observed [|T^ that the correctly defined 
lowest degenerate boundary operators introduce shifts of the boundary parameter s s±ib 
or s — > s ± ib~^ at the points of the boundary where they are inserted. Therefore, in 
order to close the conformal bootstrap, one is led to associate with a boundary operator 
two independent boundary parameters, sieft and Sj-ight, labeling the Liouville boundary 
conditions on both sides. FZZ derived a simple functional equation for the boundary two- 
point function (see the concluding section of which involves a shift of one of the 
two boundary parameters as well as the Liouville charge of the two boundary operators. 
Liouville fields. The boundary two-point function can be obtained as the unique solution 
of this equation that satisfies the refiection and duality symmetries. 

The functional equation of FZZ can be written as a difference equation with respect to 



one of the boundary parameters. In it was found that all boundary Liouville structure 
constants satisfy similar difference equations. Interestingly, these equations have a natural 
interpretation in the discrete formulation of 2D QG in terms of a gas of loops and lines 
on a randomly triangulated disc. In the discrete approach the difference equations are 
satisfied by the n-point boundary correlation functions with any n, while in pure Liouville 
theory this is so only for n < 3. Of course, this is so because starting with 4 operators, 
the contributions of the Liouville and matter sectors do not factorize. 



The general validity of this suggestion, which was proven in [19| only for the two-point 



function, follows from the formula for the three point function derived subsequently in |2C]. 
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Our original motivation for this work was to give a continuous derivation of the dif- 
ference equations for the boundary correlators obtained in ||2^,^ within the discrete ap- 
proach. This naturally led us to consider a deformation of the boundary ground ring 



structure by both bulk and boundary cosmological terms. Using the results of [19,23 



and extending the technique developed in [11], we derived a set of functional recurrence 



equations for the n-point boundary correlation functions that generalize the recurrence 
equations of Bershadsky and Kutasov [11 on one hand, and the difference equations ob- 
tained in [|19[ and ||2^,^ on the other hand. These equations completely determine all 
boundary correlation functions in 2D QG. 

In this paper we consider the realization of 2D string theory on a FZZT brane in which 
the matter field is a gaussian field with pure Neumann or Dirichlet boundary conditions. 
The target space of the gaussian field represents the Euclidean time direction x. The 
generators of the boundary ground ring, being essentially products of Liouville and matter 
boundary degenerate fields, introduce shifts both in the Liouville boundary parameter 
s and in the time direction x. The time shift is equal to the critical distance Ax = 1 
associated with the Kosterlitz-Thouless point. In general the ring relations yield functional 
equations which involve shifts both in the Liouville and matter boundary conditions, but in 
the simplest case we are considering the correlation functions depend only on the Liouville 
boundary parameter s. 

The minimal configuration of FZZT branes that closes under the action of the bound- 
ary ground ring is given by an array of FZZT branes spaced at the critical distance Ax = 1 
in the time direction. We propose a matrix model description for such a background based 
on a variant of the A^o matrix model. In this matrix model the boundary ground ring has 
explicit realization as the ring of the polynomials of of two complex matrix fields. 



2. Some preliminaries 

2.1. The boundary CFT for the FZZT brane 

We consider the non-unitary realization of Euclidean 2D string theory by a Liouville field 
(p and a gaussian matter field x with background charge cq. The background charge is 
normalized so that the conformal anomaly of the matter field is 

c=l-6eg. (2.1) 



The disc partition function of string theory with FZZT-type boundary condition [ |TP| , |23| is 
defined by by the following effective action, which comprises a bulk and a boundary term: 



A[xA] = I (J^[{V<pr + iVx? + {QcP-ieox)R]+^e''^^ 

M 

+ j - ieQx)K + iiB e'"^^ + ghosts 



(2.2) 
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where the couphngs ji and jiB are are referred to as bulk and the boundary cosmologi- 
cal constants. The background charges are expressed in terms of the Liouville couphng 
constant h as 

Q=T+h, eo = --b. (2.3) 



With the choice (2.3) the two background charges satisfy Q'^ — Cq ~ 4, which is equivalent 
to the balance of the central charge ctot = + + Cghosts = (1 + 6Q^) + (1 — 6eo) — 26 = 0. 
We will consider the generic situation when b is not a rational number. 

After mapping the disc A4 to the upper half plane {Imx > 0}, the curvature term 
disappears 

oo 

^[0,X]= J d^x(J^[{Vcl>f + {Vxf]+fie^''l'^ + J dxfiBe'^ + ghosts (2.4) 

lma;>0 — oo 

and the background charges are introduced through the asymptotics of the fields at spatial 
infinity 

0(x, x) ~ — Q log x(x, x) ~ — Co log (2.5) 

The boundary term encodes a inhomogeneous Neumann boundary condition for the Liou- 
ville field 

i{d - d)4>{x, x) = 47r//B e^-^^^'^^ {x = x), (2.6) 

which describes the FZZT brane characterized by the the parameter /ib- The matter field 
is assumed to satisfy pure Neumann boundary condition: 

i{d-d)xix,x) = (x^x), (2.7) 

which is of course Dirichlet boundary condition for the T-dual field x{x,x). In this paper 
we consider the simplest situation where no matter screening charges are added. 

2.2. Bulk and boundary vertex operators (closed and open string tachyons) 

The bulk primary fields V^^\x.iX)^ or left/right moving on- mass-shell closed string 
tachyons, are defined as the bulk vertex operators 

V(+) = -7(6P) e»(eo-P)x+(Q-^')0 

T (2.8) 

P ^ \ h } 

where 7(0;) = r(a;)/r(l — x). This normalization removes the external 'leg pole' factors in 
the correlation functions and is the one to be used when comparing with the microscopic 
realizations of 2D QG. 
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We will be mainly interested in the correlation functions of the boundary fields that can 
be inserted along the boundary of the world sheet. The boundary primary fields B^^\x), 
or left/right moving on- mass-shell open string tachyons, are defined as the boundary vertex 
operators 

= 1 T{2bP) e^(^-o-P)x+(^Q-P)<^ 

T (2.9) 
gi-) ^ i p (-fP) e'(^^°"^)^+(^'^+^)'^. 

As any CFT boundary operator, the open string tachyon is unambiguously defined only 
after both left and right boundary conditions are specified, which should be done both for 
matter and Liouville components. 

In a theory in which the ends of the open strings freely propagate in the two- 
dimensional space-time (/is = 0), both fields satisfy pure Neumann boundary conditions. 
In a theory with FZZT branes, the Liouville left and right boundary conditions are de- 
termined by the values si = s^'^^^ and S2 = s"^^^^^ of the uniformization parameter defined 



in (1.1). Slightly modifying the notations of FZZ [19|, will denote such an operator as 
^'^[B^^^Y^. Geometrically the boundary fields create open string states whose ends propa- 
gate in two different FZZT branes labeled by si and S2- 

The physical observables can be represented in two pictures: either as (1,1) -forms and 
integrated over the the world-sheet (for the bulk fields) and 1-forms integrated over the 
boundary (for the boundary fields), or as BRST-closed 0-forms: 

/ (fxV^J^\x,x) ccVl^\x,x) 

^ (2.10) 

dxB^^\x) ^ cB^^\x). 

Here and in the following b, c and b, c are the reparametrization ghosts respectively in the 
bulk and on the boundary. The second representation is more appropriate as far as the 
ground ring structure is concerned. 

2. 3. Reflection property 

The tachyons of opposite chiralities are related by the Liouville bulk and boundary reflec- 
tion amplitudes (see Appendix A) 

(2 11) 

^-[B^^Y=D^P~\suS2) ""^[B^p Y 
2.4- Physical states in 2D QG 

For each value of the momentum there is only one vertex operator that corresponds 
to a physical state. The physical operators are 

Vp=^^fl (^>0), Bp^l^n (^>0) (2.12) 

^ vjr^ < 0) 1 B^p^ (p < 0). 
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The "wrongly dressed" operators are related to the physical ones by the Liouville reflection 
amplitude. 

2. 5. Degenerate fields 

By degenerate flelds in Liouville quantum gravity we will understand the gravitationally 
dressed degenerate matter flelds. In our case these are the on-shell vertex operators (2.12) 
with degenerate matter components 

Vrs=Vr/b-sb. S.s = (r,sGlN). (2.13) 

The fusion rules for these flelds are the same as the fusion rules for the degenerate flelds in 



the matter CFT . Note that the Liouville components of these flelds are not degenerate 
Liouville flelds. 

The operators that span the ground ring are degenerate with respect to the full con- 
formal algebra. Such an operator is constructed by applying a rising Virasoro operator of 
level rs — 1 to a product of matter and Liouville degenerate flelds with Kac labels r, s. The 
conformal weight of the product is zero: 

2. 6. Normalization of the hulk and boundary cosmological constants 

It is also convenient to redeflne the cosmological constants // and /is according to the 
normalizations (2.8) and (2.9) of the vertex operators. The new bulk and boundary cos- 
mological constants A and z are deflned as the couplings of the operators vi^^ and Bl~^j2 
respectively: 

f^e''^ = AVi+\ e'^ = z b[^]^. (2.15) 

This gives 

-=7r7 62, — = 2.16 

(1 hb r(i-6^) 

The uniformization map (1.1) now reads 

2 = Mcosh(7r6s), M = VX. (2.17) 

Here we introduced the constant M whose physical meaning is that it is equal to minus 
the boundary entropy due to the fluctuations of the Liouville flelJ. 



^ When comparing with the discrete approach, the cosmological constants A and z are to be 
identified with the corresponding quantities in the microscopic formulation of the strings defined 
on Dynkin graphs [27|. The constant M equals half the loop tension in the loop gas description 
of 2D string theory [27|,p8|. The uniformization parameter s is defined as in [p^; it is related to 



the parameter r of |27] by r = irbs. 
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2. 7. The self-duality property of Liouville quantum gravity 

All correlation function in Liouville theory are invariant w.r.t. the substitution 



19 



b^b=l/b, A^A, M^M 

where 

A = A^/''', M = M^/^\ (2.18) 

The boundary parameter s is self-dual. The duality transformation for the correlation 
functions follows from the one of left and right chiral tachyons 

Vj,^)=v5.\ '^[B'^^^Y' = ''[B^J)Y\ (2.19) 

where the bar means complex conjugation. 

It is clear from (2.19) that the non-unitary CFT under consideration is not self-dual 
in the strict sense. The duality transformation interchanges the states of a pair of CFT 
characterized by matter background charges cq = ^ — b and cq = b — ^ = —cq. The 
duality becomes a symmetry (a chirality flip) only when it is accompanied by a complex 
conjugation. 



3. The (bulk) ground ring 

We start with a short review of the ground ring structure for a theory defined on the sphere 
{dAi = 0). The ground ring operators are obtained by applying raising operators of level 
rs — 1 to the product of two degenerate matter and Liouville fields with Kac labels r, s. 
The resulting operators have conformal weights A = A = 0. The ring is generated by the 
lowest o two perators [H 

a+ = - |bc - bd,{<P - tx)\^e-'''^^+'^^ 

a- =-\hc-b-^^,{(t> + ix)\^e-^^^-'^'^■ 
where b, c are the reparametrization ghost and anti-ghost fields. The ground ring is 
spanned on the polynomials (a-|-)"^(a_)'^ with m, n G In the case of non-rational 

6^ the ground ring contains no other relations and has an infinite number of elements 
labeled by the integers r, s > 1. 

A crucial property of the operators a± is that their derivatives dj:a± and d^a"^ are 
BRST exact: d^^a- = {Qbrst, b_ia_}. Therefore, any amplitude that involves a± and 
other BRST invariant operators does not depend on the position of a±. This property 
allows to write recurrence equations for the correlation functions from the OPE of a± and 
the other BRST-invariant fields. 
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3.1. The limit of of free fields (K = Q) 

The vertex operators (2.8) form a module under the ground ring: 

a+vi,+) = -V^+_),, a.V^p-^ = -V^p-\. (3.2) 

and also 

a+vjr^ = a_v{,+^ = 0. (3.3) 

Both relations (3.2) and (3.3) follow from the free field OPE§ and are true up to commu- 
tators with the BRST charge. While the first one survives in the correlation functions, 
the second one receives non-linear corrections. Due to the contact terms, the last relation 
should be modified in presence of an integrated vertex operator to 

a+V^T^ j (fz V'p-^ = V'flp^^, (3.4) 



These relations imply a set of recurrence equations |P,^,|TT|, which determine completely 
the resonant tachyon amplitudes (see Appendix B). 

3.2. A ^ 

The first relation (3.2) survives the perturbation with A 7^ 0, while the second relation 
(3.3) gets deformed. The deformation is linear in A and can be calculated pertubatively 
by considering a first order insertion of the Liouville interaction — /i J e^^*^ = —A J Ve^'* 
and then applying (3.5) : 

a_vj,+^ = -A V},+_V (3.6) 

This procedure has been justified by various self-consistency checks in Liouville theory 
|]T^,|TB[. By the same argument one finds the action of the operator a+. This time one 
should consider a linear insertion of the dual Liouville interaction —jl J e^'^^^ = —A f V^J^y 
The result is 

a+vlr^ = -A V^/„ A = A'/'\ (3.7) 

The deformed ring relations allow to derive an (over determined) set of recurrence equations 
for the bulk correlation functions (see Appendix B). 

Here we consider the simplest situation where the matter field obeys U{1) fusion rules. 
In general the free-field representation of 2D QG should be completed by adding to the 
effective action (2.2) the two matter screening charges 



-2ibx 



The normalization of the fields is determined by the action (2.4) and is the same as in [19|. 
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Then eqns. (3.6) and (3.7) will further deform to 



a-V^^^ = -A V^^\ + V^^j, 



a+V^-^ = -A V^/, + V;-_\/,. (3.9) 



The ground ring structure in presence of screening charges has been considered in |3C,31 . 



3.3. The case b=l 

The case 6 = 1 plays a special role [^. Consider the action of the product a+a_: 

When 6 = 1 we get, since A = A and cq = 0, 

a+a_=A (6=1). (3.10) 

This equation has a direct interpretation in the profile of the Fermi surface in the cor- 
responding large matrix model 0. This analogy can be pursued further to a generic 
perturbation by bulk tachyons and relate the ground ring structure and the Toda integrable 
structure in the matrix realization of the c = 1 string theory |32-37[ . 



4. The boundary ground ring 

4.1. The boundary ground ring for pure Neumann boundary conditions (hb = 0) 

We first introduce the boundary ground ring for 2D open strings with of pure Neumann 
boundary conditions for both coordinate fields, (p and following Bershadsky and Kutasov 
[Til] . The two generators are defined, following the same logic as in the bulk 



A+ = - [be - lbd^{<f> - ix)] e-^'~'^^+^^^ 



2 (4.1) 
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A_ = -[be - -6 ^d^{(i) + ix)] e ^ 

The two operators are related by a duality transformation combined with complex conju- 
gation: 

i+ = A_, i_=A+. (4.2) 
9 



They are BRST closed: dj:A± = {Qbrst, b_iA-|-} and have A = 0. The open string 
tachyons (2.9) form a module with respect to the ring generated by these two operators. 
The action of the ring on the tachyon modules is generated by the relations 

= B^+l^, A_B^-^ = B^-\ (4.3) 

and 

A+B^~^ = A_b'^^^ = 0. (4.4) 

If one exchanges the places of A and V, there will be minus on the rhs, since A is odd and 
V is even wrt x —>■ —x: 

B^+^A+ = -B^+l^, B^-^A. = ~B^-\. (4.5) 
Again, the first relation (4.3) is exact and the second relation (4.4) gets deformed in 



presence of integrated boundary tachyon fields 



A_B^^^ I dxB^+^ = B^^l^ , . (4.6) 

Since on the boundary the ordering is important, there is a second relation 

[dxB^;^A.B^;^= ^'pIp - (4-7) 

The coefficients in (4.6) and (4.7) are obtained as standard Coulomb integrals (see for 
example |]l^). Similarly, one finds for A^ 



and 

,(-), .(-) sinf(Pi + P2) ,(-) 



dxB''A+B'' = ^^^^ :y^B''^^,. (4.9) 

+ sin^P2 sin^Pi Pi+P2+^ ^ ^ 

Using these relations, Bershadsky and Kutasov obtained recurrence equations for the n- 
point open string amplitudes (Appendix C). 

4-2. The boundary ground ring in presence of FZZT branes 

The ring relations can be generalized to the case of nonzero boundary cosmological constant 
using the results of the paper by V. Fateev, A. Zamolodchikov and Al. Zamolodchikov 
19(1 . The crucial observation made by the authors of |T^ and later confirmed by B. 
Ponsot and J. Teschner |2C] is that a level-n degenerate boundary Liouville field e""^"^ 



has vanishing null-vector and therefore a truncated OPE with the other primary fields if 
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either sieft — bright = ibk or sieft + Sright = ibk, with k = — n/2, — n/2 + l,...,n/2. By 
the duahty property of Liouville theory, the degenerate boundary fields e~'^'^^^ exhibit a 
similar property with h replaced by 1/6. No direct proof is supplied for this statement, 
but it was shown to be consistent with the exact results obtained in boundary Liouville 
theory. In particular, the above property leads to a pair of functional equations for the 
Liouville boundary refiection amplitude, whose unique solution coincides with the result 
of the standard conformal bootstrap. 

According to [M, the operators (4.1) should be defined byS 



A 

and the relations (4.3), (4.6), (4.8) and (4.9) understood as 

''[B^~^Y[A_Y^'^ = - '^[B^~},Y^''' (4.11) 

^ 2 

etc. As in the bulk case, the relation (4.4) is deformed by a linear insertion of the boundary 
Liouville interaction. However, in this case the situation is more subtle since the boundary 
interaction depends on the point it is inserted. We use the fusion rules that follow from (4.6) 
and (4.7) where one of the operators is the Liouville boundary interaction —z{s) x J B^^j^. 
The result is 

sij^_|si±i6jg(+)|S2 = c±x '^[B^+},Y^ (4.12) 

^ 2 

with the coefficient given by the sum of the contributions of the three possible insertions 
of —z{s) X J i3e^'* with s = Si, Si ± ib and S2 (fig-1): 

C, = -M (cosh(....) .7 f + ^' -^ '"-"'-^t;' ^ '^'i + ^g!^) . (4.13) 

\ sm(7roeoj sm(27roi^J sm(7roeoj sm(27r0-rj / 

The coefficient C± can be nicely written as 

cosh[7r6(si ± 2zP)] + cosh(7r6s2) 

C± = —M — — 

sm27roP 

cosh [nb (^4^ ± iP)] cosh [nb ± iP)] ^^'^^^ 



-M- 



\ sin 2'KbP 



^ We omit the boundary conditions for the matter field. In our case of gaussian field x with 
pure Neumann boundary conditions, those are given by the values of the dual field x on both 
sides. The momentum conservation yields Xright — Xieft = i??^^. In calculation expectation values, 
the final result will not depend on the only residual parameter, the global mode of x- However, if 
screening charges are allowed, the matter boundary conditions matter. 
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Fig. 1: The three possible insertions of the boundary Liouville interaction. 



Similarly one finds for the action of on Bp 



(-) 



(4.15) 



with 



~ _ cosh[7r (gj ± 2zP) /6] + cosh(7rg2/^) 

^~ sin2fP 

^ cosh [f ± zttP)] cosh [f ± zttP)] 

isin(2fP) 



(4.16) 



Note that the fusion coefficients are symmetric with respect to the change the sign of any of 
the boundary parameters. This is a general property of all correlation functions and comes 
from the fact that the s-space is an orbifold: the points s and —s should be identified. 



4-3. Functional equations for the correlation functions of boundary primary fields 
Consider a boundary correlation function of the form (fig. 2) 



Wp,,p,Usi, s, S2, S3, ...) = ( [Bp^ViBp^r [Bp,Y\.. ) (4.17) 

The realization of the physical boundary fields depends on the sign of the momenta and is 
given by (2.12). The amplitude (4.17) is by construction symmetric in its arguments. It 
is analytic in the boundary parameters s, si, ... but not in the target space momenta. By 
the momentum conservation the correlation function is zero unless X]a;(i^o ~ -Pfc) = ^o- 
Equation (2.19) implies the identity 

Wp,,P,,... = W_p„_p„... = Vr_p„_p„.... (4.18) 
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Let us assume for definiteness that Pi < and P2, -P3 > 0. Then the amphtude (4.17) 
is reahzed as 

i^p„p„P3,...(«i, «2, «3, ...) = ( [4;^]^[4^ [B^pjv'- ) (4-19) 

By the symmetry (4.18) the amphtude can be written as 

Wp„p„P3,...(si, .2, S3, ...) = ( ...^n^U]^^[^yj1^i"'A]^^ ) • (4-20) 
We will use these two representations to derive two independent functional identities for 




Fig. 2: Our convention for the ordering of the boundary parameters and momenta. 

Consider the auxiliary correlation function F with an operator A- inserted in the r.h.s. of 
(4.17): 

Let us assume that Pi < — | and P2 > |. Then, evaluating F by (4.11) and by (4.12) and 
equating the resultsl we get the following functional identities for the correlation functions 
of three or more fields: 

sin (27r6P2) VrPi,P2,P3,...('Si, « ± S2, S3, ...) = 

2Mcosh [nb ± zPs)] cosh [nb ± ^^2)] W^p, + |,p2-|,p3,...(si, s, S2, S2, ...) 

+ M^P, + |,P,+P3-^,...(S1, s, S3, -)■ 

(4.21) 

The last term represents a correlation function with one operator less. 

A dual equation can be obtained in the same way by evaluating the complex conjugated 
function F using the representation (4.20) and the relations (4.2): 

sin (^^2) W^Pi,P2,P3,...(si, s ± f , S2, S3, ...) = 

2MV^^ cosh [f ± ZP2)] cosh [f ± ZP2)] Vrp^ + ^,P,_^,P3,...(S1, S, 82, S2, ...) 

+ W^p, + ^,p,+P3_|,...(si, S, S3, ...), 

(4.22) 

^ Our notations do not distinguish between integrated and non-integrated fields; the rules are 
the same as in Appendix C. 
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where we assumed that -Pi < — ^ and P27-P3 > These identities form a set of (over- 
determined) functional recurrence equations that generahze the recurrence equations of 

Taking the differences in (4.21) and (4.22) one finds a pair of homogeneous difference 
equations which do not have analog in the = limit. We write them in operator form: 



sm{bds) -Msinh(7r6s) 65(^^1-^^-2)1 VTp^ (s^, s, 53, S3, ...) = 



(4.23) 



and 



sin(ic>,) - M^/'^' sinh(f s) e^^^^i-^^^)] Wp,,p,,p,,..Xsu s, S2, S3, •••) = 0. (4.24) 



where it is assumed that Pi < — | and P2 > |- The difference equations (4.23) and (4.24) 
have been first obtained in the microscopic formulation of 2D QG as a gas of loops and 



lines on the randomly triangulated disc in p^j2^ . 

Functional equations for all values of the momenta can be obtained by the following 



rule ||2J,^. If the sign of one or both momenta changes after the shift, one should apply 
the refiection property with the amplitude -D^ •* (si, S2) given in Appendix A. The details 



are explained in [p3 



The functional equations (4.21) and (4.22) can be used to evaluate the two- and one- 
point functions as well as the disc partition function on the FZZT brane. These quantities 
are difficult to calculate directly in CFT because of the subtleties associated with the 
residual global conformal symmetry. The calculations are presented in Appendix D. 



5. Matrix model description 

The old idea of Polchinski that matrix models for non-critical strings describe backgrounds 
populated by a large number of unstable D-branes was recently given a concrete shape in 



The boundary state of such a D-brane is a product of the ZZ boundary state 
for the Liouville field, localized at large (p, and a boundary state for the matter field 
corresponding to some local boundary condition. For example the standard matrix model 
for the c = 1 string, the Matrix Quantum Mechanics (MQM), describes the case of pure 
Neumann boundary condition for the matter field (or equivalently pure Dirichlet boundary 
condition for the dual field). 



In this section we will show, using the results of |]2^,^, that there exists a neat 



matrix model interpretation of the boundary ground ring relations. The matrix model 
dual to the world sheet theory (2.2) is an infinite matrix chain, which can be viewed as 
a special discretization of MQM. The logic is the following. First, as it was shown in 
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|2^,^, the difference equations (4.23) and (4.24) have straightforward interpretation in 
the loop gas formulation of non-critical string theories [ ^p7| , |28[1 in which formulation the 
string path integral is defined microscopically as a sum over planar graphs embedded in a 
Dynkin diagram of ADE or ADE type. This allows to find the target-space equivalent of 
a degenerate boundary field. Second, the loop gas formulation can be restated in terms of 
the ADE and ADE matrix chains studied in p9| - |5^ . The matrix model description of the 
world sheet theory (2.2) is provided by the A^o matrix chain. 

It is natural to think of the Aoo matrix chain as a marginal perturbation of MQM 
by a periodic potential ~ cos(27rx) with the effect that each DO-brane decomposes into 
an array of D-instantons associated with the minima a; G Z of the potentiafl. The initial 
state of such a system is an array of FZZT branes placed at the sites of the A^o Dynkin 
graph (A^ branes per site). 



5.1. Loop gas and branes 

The continuum limit of the loop gas model can be given the following interpretation within 
the 2D string theory. It describes the evolution of closed strings in presence of an array of 
FZZT branes spaced at the critical distance Sx = 1 in the Euclidean time direction. The 
background charge for the matter field is introduced by boosting the array of D-branes 
with momentum cq. The target space of such a string theory is a direct product of the 
Liouville direction and the discretized time direction 



X — nxb, X G Z, 



(5.1) 



in which only the points of the array are retained. The string world sheet is divided into 
domains associated with different branes. It is assumed that only jumps a; — x±l between 
two nearest- neighbor branes are allowed. This leads to the geometrical picture of a gas of 
non-intersecting loops, representing the domain walls on the world sheet (fig. 3). 




Fig. 3: A disc world sheet with boundary in the x-th brane. 
Another possible matrix model description of such a perturbation is proposed in |55]. 
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Creation of a domain wall costs energy which is proportional to the length of the domain 
wall and which is taken into account by the loop tension M. The statistical weight of each 
domain wall is e~^^^ where i = f e^"^ is the length of the domain wall. The loop gas 

domain 
wall 

describes a sector of the world sheet theory (2.2) containing certain class of degenerate fields 
(2.13) (for a recent discussion see [0). Let us remind that the field x satisfies Neumann 
boundary condition and that the dual field x satisfies Dirichlet boundary condition and 
that the lowest degenerate boundary operator B12 is realized as a magnetic operator with 
magnetic charge P = 6/2. The operator Bi^2 creates a jump of the boundary to an adjacent 
brane in the x space. Geometrically this operator is described by a line starting at some 
point on the boundary. Similarly the degenerate boundary operator Bi^n+i with magnetic 
charge P = nb/2 creates a jump by n steps and is described by n lines starting at given 
point on the boundary. The correlation function of two such operators is a sum of world 
sheets as the one shown in fig. 4. 




The degenerate boundary operators Bn+1,1 with n = 1,2,... can be constructed as 
electric operators interpolating between two Neumann boundary conditions for the field 
X- The Neumann boundary condition is characterized by the global mode representing 
the increment X — J dx along the boundary. The action of these operators closes in the 
discrete set of boundary conditions characterized by 

X = TTx/b, X G (5.2) 

The duality symmetry b ^ 1/b exchanges the x and x spaces. In the dual loop gas 
model the time coordinate x is discretized with spacing Ti/b. The weight of a domain wall 
is e~^^^^, where i = f e^/'^dx. The dual loop tension M = M^/^ is the dual boundary 

^ domain 
wall 

cosmological constant for the theory (2.2). 



16 



Now let us see what is the meaning of the operators of the boundary ground ring 
in the loop gas model. Geometrically the degenerate boundary operator with momentum 
Bn+i,i is obtained from Bn,i by adding an extra line. This reminds the fusion rules for 
the matter fields. However, the process of adding a line cannot be interpreted directly as 
a fusion of two boundary tachyon operators because the Liouville dressing factors obey a 
different fusion rule. The correct interpretation comes from the fusion rules (4.3) involving 
one tachyon and one operator of the boundary ground ring. Thus the operators and 
A_ generating the boundary ground ring have the geometrical meaning of the operations 
of adding or removing lines. The action of these operators is encoded in the left and right 
boundary conditions for the matter field: 

A+ = '^'""{hA+W A_ = 

(5.3) 

A+ = '^^'/'ra+r, A_ = ^-t'tP-l-. 



5.2. A self-dual Aoo matrix chain 



The matrix chains associated with Dynkin diagrams were originally defined in terms of 
a set of N X N hermitian matrices associated with the nodes x of the Dynkin diagram 
and N X N complex matrices A(^x,x') associated with the links {x,x'). All these matrix 
chains have a loop gas description and therefore a direct Coulomb gas interpretation. The 
complex matrices A(^x,x') represent the one type of degenerate boundary fields (only electric 
or only magnetic ones). 

Each value of the central charge is realized by two such matrix models, related by 
duality. For example the theory with b = and central charge c = 1 — p(^p^i-^ is realized 
either as an Ap chain in the dilute phase or as an Apj^i chain in the dense phase. The 
electric degenerate boundary fields have a good description in the Ap matrix chain while 
the magnetic degenerate boundary fields have a good description in the Ap+i matrix chain. 

Remarkably, each such pair of matrix chains related by duality can be formulated 
as a single self-dual chain of complex matrices by the following simple prescription. The 
hermitian matrix variables associated with the nodes of the Dynkin diagram should 
be considered as composite fields B^bI or bIb^, where B^ are N x N complex matrices. 
Here we will consider only the Aoo chane, which provides the matrix model description of 
the theory (2.2). The case of the rational ADE theories will be studied elsewhere . 



The partition function of the self-dual A^o matrix chain is defined in terms of the 
N X N complex matrices A(^x,x') = and B^ (x G 2): 

2= [U dAMdBxdBl e-^(^'^"^'^'^"^) (5.4) 
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where the simplest action with linear potential for ylylt and IS 

S{A, A^B, st) = J] It (2TA,aI + 2TB,bI - B.bU^aI - 4+i5x+i4^x) • 

(5.5) 

(Of course this integral has only formal meaning because of the infinite volume of the 
chain. It can be regularized as the limit n — > oo of the periodic chain An-) 

Alternatively one can write the partition function as an integral over the positive 
definite hermitian matrices 

^^ = B^bI ^^ = A^aI (5.6) 

and the unitary matrices Uj;. Let us denote by [rf$]>o the flat measure over the positive 
definite hermitian matrices and by [dU]^^^^^^ the Haar measure on the group U{N). Then 
the integral (5.4) can be written as 

2= [U [d^.Uid^.Um^U,^, e-^(*'*'^) (5.7) 

The integration over the unitary group can be done explicitly and the partition function 
can be written as an integral with respect to the eigenvalues (pxj and (pxj of $c and 
respectively. Using the Harish Chandra-Itzykon-Zuber formula we write (5.7) in the form 

oo 

2^ Jjl dcpxj d^xj e-^^'^^-^-'^^^-'^ dete'^-'^^^'" dete'^'^'^^^+^'K (5.9) 

We are interested in the 't Hooft limit N ^ oo with k = N/T"^ fixed. Then the 
ground state of the matrix model is characterized by the classical spectral densities Px{(l>) 
and Px{(t>) of the hermitian matrices (5.6). The 't Hooft limit exists for k < Kc, where Kc 
is the critical coupling where the planar graph expansion diverges. One can show that in 
this case Pxi^P) and Pxi^P) vanish for (pxT^Px ^ T, up to exponentially small in terms. 

We will be interested only in the scaling limit, which describes the infinitesimal vicinity 
of the critical point. It is technically simpler to treat the parameters N and T as large but 
finite. Then for given T the critical point is achieved at Nc — KcT"^. The value Nc plays 
the role of a cutoff parameter, which we eventually sent at infinity, while the cosmological 
constant is p, ^ Nc — N. 

The matrix chain (5.5) describes the particular case 6 = 1 of the world sheet theory 
(2.2), when the theory is unitary. A non-zero background charge eo = 1/b — b can be 
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introduced by imposing by hand a a;-dependent background. Then the spectral densities 
are solutions of the saddle point equations under the constraints 

pM = e'""''^/' p(<^), pA<P) = p{<l>). (5.10) 



o 



Relation to the standard A^o chain ^4M1 



One can formally integrate with respect to the (^-variables in (5.9) neglecting the contri- 
butions of infinity. In this way one neglects exponentially small in terms, which do not 
influence the genus expansion. The result of the integration is the partition function of 



the hermitian A^o matrix chain [49 



^'-^ 



By integrating first with respect to the ^-variables one obtains the same eigenvalue integral 
for the dual field ^x'- 



^ VVj.k ~ Vx,j - (Px+l,k^ 



^'-^ 

The saddle point equations for these integrals lead to spectral densities Pxifp) ^-nd Px{(t>) 
that vanish for > T, which justifies the fact that we neglected the contributions at 
infinity. 

o Relation to Matrix Quantum Mechanics 

It is known [p7[ that the hermitian Aoo matrix chain (5.11) and MQM share the same 



x-space correlation functions of on-shell closed string tachyons. Thus the Aoo chain is 
expected to yield a discretization of the singlet sector of MQM at euclidean time distance 
P = n/2 in the normalization in which the self-dual distance is /3 = 27r. Let us see that 
this is indeed the case. 

The singlet sector of MQM is equivalent to a system of nonrelativistic fermions in 
upside-down oscillator potential U{y) = —^y^, representing the eigenvalues of a hermi- 
tian matrix variable (see the review 0] and the references to the original papers therein). 
The evolution of each eigenvalue y in MQM is described by a one-particle Hamiltonian 
Hq = — |(9y + y"^). The heat kernel for this Hamiltonian is given by 

Kf3{y,y) = y e "\y) = . exp — -— . 5.13 

V27rsmp \ zsmp / 
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Therefore the Laplace kernels entering the determinants in (5.9) can be viewed, after a 
linear change of variables 

y = 0-T, y = 0-T, 

as evolution kernels (5.13) with {3 = ^tt. Let us denote by 7i = Hq the MQM Hamil- 

tonian in the singlet representation and by Vt the projector restricting the integration to 
the semi-infinite interval [—T, oo] : 

N 

er = l[e{yi-T). 

Then the partition function (5.9) for a periodic chain An can be written as the trace 

e-^'^^j . (5.14) 

The projector has the effect of an infinite potential wall placed at distance T from the 
origin, which stabilizes the inverse oscillator potential. For given T the Fermi level grows 
with the number of the eigenvalues N and at some critical N = Nc reaches the top of the 
potential. The cosmological constant /j, Nc — N measures the deviation from the critical 
point. 

5.3. Disc loop amplitudes in the matrix chain 

The classical background of the matrix chain is determined by the (unnormalized) disc 
loop amplitudes 

W,(^) = /Tt 

\ z — BxB. 

and 

t^,(5) = /tt- 1 ^ 

The loop amplitude (z) represents a meromorphic function in the spectral parame- 
ter z with a cut along the support [0, a] of the eigenvalue distribution, where < a < T. It 
is determined from the saddle-point equation for the integral (5.11), which can be written 
in the form of a boundary condition along the cut: 

W^{z + iO) + W^{z-iO) + W^+i{2T-z) + Wa,-i{2T-z) = 2T {0<z<a). (5.17) 

The dual loop amplitude Wx (z) has the same properties with z, a replaced by z, a. 



d4>px{4>) 



z — 



d(t)px{(t') 

Z — (j) 



(5.15) 



(5.16) 
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The true ground state does not depend on x and describes the unitary theory with 
matter central charge c = 1. The non- unitary theory with background momentum cq is 
obtained by imposing by hands the foUowing constraints on the loop amplitudes:! 

W^{z) = e'^o^/'' W{z), W^{z) = e-^^°^^ W{z). (5.18) 

Then the equation (5.17) becomes 

W{z + iO) + W{z -i0)-2 cos{ti/}?) W{2T - z) = 2T {0 < z < a). (5.19) 

This equation is readily solved in the scaling limit, when "^^"^'^ and -^^^^ are considered 
small. Then the branch point at z — T can be resolved by the parametrization 

z-T = M cosh(7rs6), M = T - a, (5.20) 

and (5.19) becomes a difference equation 

W{s + ib) + W{s -ib)-2 cosh n/b'^ W{s) - 2T. (5.21) 

Prom here we find 

W{s) = ciT + C2T^-^^ M^/^' cosh(7rs/6). (5.22) 

where ci and C2 are 6-dependent constants. The solution reproduces, up to a linear trans- 
formation, the loop amplitude (D.16) in the world sheet CFT. 

In the following we will retain only the scaling parts of z and W, assuming a normal- 
ization that matches with the result from the continuum CFT: 

z = Mcosh(7r6s), W= -— cosh(7rs/6). (5.23) 

^cos(^) 

In this normalization the eigenvalue density is given by 

p{z) = ' + = bM'/"' sinh(vr./6). (5.24) 

Proceeding in the same way we find the scaling part of the dual loop amplitude (5.16) 

z = Mcosh(nb/s), W = -— j- cosh(7rs6), (5.25) 

ocos(7ro^) 

where M = T — a. 



^ It is actually sufficient to impose only the second condition; then the first one follows from 
the correspondence W z, z ^ W . 
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The loop amplitudes (5.23) and (5.25) are related by a duality transformation 

Vr = - , z = -lsmi^)W. (5.26) 

sm nb^ ° " 

Intuitively, the duality in Liouville gravity is a kind of electric-magnetic duality on the 
world sheet, whose target-space equivalent is a Fourier transformation for the eigenvalues 
of the random matrix. 

In the planar limit the matrix model is described by a classical Hamiltonian system 
in which the spectral parameter z and the resolvent W form a pair canonically conjugated 
variables: the coordinate and the momentum. Then the duality make sense of a canonical 
transformation exchanging the momentum and the coordinate. The meromorphic functions 
W{z) and W{z) give two parametrizations of the complexified classical trajectories of the 
Hamiltonian system. 

The physical sheets of the Riemann surfaces of W{z) and W{z) can be viewed, together 
with the identification (5.26), as two charts defining a complex curve with the topology 
of a sphere with two punctures. In the rational models of quantum gravity the Riemann 
surface is given by an algebraic curve. In the case of generic b^ the equation of the algebraic 
curve can be written as a functional equation 



W'^iz) + W^{-z) - 2 cos(^) W{z)W{-z) = b'^M^/^^ (5.27) 



or its dual 



W^iz) + W^i-z) - 2cos(7r62) W{z)W{-z) = ^ M'^^ . (5.28) 



These equations are particular case of more general functional equations that hold before 
the scaling limit and allow to express the positions a and a of the branch points in terms 
of T and A^. In particular, the non-scaled functional equation for W{z), whose derivation 
presented in Appendix E, implies the scaling M'^^ ~ A. 

5.4- Planar graph expansion and mapping to the world-sheet CFT 

The perturbative expansion of any observable in the theory with action (5.5) is a sum of 
planar graphs composed from the vertices 

A. \ aj r 



x + 1 



where the continuous and dashed oriented lines represent respectively the propagators of 
the A and S-matrices. The planar graphs can be considered as discretized world sheets 
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embedded in TL. The loops made of the A^j-propagators can be viewed as domain waUs 
separating domains with time coordinates x and x + 1. 

We are interested in the disc amphtudes representing various boundary correlation 
functions of tachyon operators. The simplest ones are expectation values of the resolvents 
(5.15) and (5.16), which are the equal to the boundary one-point functions in the world 
sheet CFT with Neumann and Dirichlet boundary conditions!. 

In the planar limit W^^ is a sum of planar graphs with the topology of a disc, made by 
A and B type lines. The A-lines can only form loops while the S-lines can either close in 
a loop or have its both ends among the external legs. The condition (5.18) is equivalent 
to associating x-dependent phase factors with the vertices and faces of the planar graph 
as explained in section 2 of [^. For example, the graph shown in fig. 5 gives a particular 
discretization of the world sheet with three domain walls shown in fig. 4. 





Fig. 5 a,b: Planar graphs contributing to the resolvents of BxBt and AxAt- 



The dual amplitude VF^ is a sum of planar graphs such that the S-lines are always 
closed and the A-lines either close or connect two external legs. In this way the complex 
matrix chain allows to define both the Dirichlet and Neumann boundary conditions for 
the loop gas in the sense of and [^. Exchanging the A and B lines is equivalent to 
the duality transformation, which exchanges Dirichlet and Neumann boundary conditions. 
In the continuum limit the two amplitudes are given, up to rescaling and subtraction of a 
constant, by the expressions (D.16) and (D.17). 

The boundary ground ring is represented in the matrix model as the algebra of the 
polynomials of the A and B matrix variables. For the generators of the algebra the 
correspondence is 

A-\- > Ax A— > Bx 

_ _ (5.29) 



^ Strictly speaking, the Neuman boundary condition is generated by the resolvent of the op- 
erator {A -\- A^)'^ . However it has the same scaling limit as the resolvent of A^ . 
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The boundary segments with boundary condition (x, s) or (x, s) are represented by 
the resolvents of the fields = B^B^ and = AxA^ correspondingly: 



z(s) - ^x 

^ \ (5.30) 

The degenerate boundary operators Bn^i and i3i „, are created by inserting in the trace 
polynomials of A and B matrices according to the following dictionary: 

X+n7r/& [^n.l] X ^ Ax-\-n---Ax x+'^t'"'' [^1 !"■] X ^ Bx-\-n---Bx 

(5.31) 

ft ft 

X-n7r/&[^-n,l]x ^ A^_^...Ax [^1, -n] x B^_^...Bx- 

For example, the boundary two-point function Dp{s, s') for P = |eo + ^ corresponds 
to the following expectation value in the matrix model (Fig. 4): 

= (x+nf [^ri+l,l]x['^-"-l -l]x+nt 

t t 1 1 \ (5-32) 

This correlation function was calculated from the matrix model side using the loop gas 
technique [p^|23[] and the result is in agreement with the prediction of the Liouville theory. 
Another example is the amplitude of two excited twist operators with momenta ±P^^ i q = 
±(n + |)/26, which interpolate between Dirichlet and Neumann boundary conditions: 



1 (n+l)('S, S') = (I \B n + l/2 Y\B r. + l/2 ]\ , ^. 

' 2/ \ A 2b ^ 2b ^ / disc 



t I 1 1 \ (5-33) 



This correlation functions has been calculated (for the case n = 0) in and given a 



Liouville CFT interpretation in |24,23 



6. Concluding remarks 

1. The ring relations and the functional equations derived in this paper can be generalized 
to the case of any matter CFT with c < 1. This can be achieved by replacing the U{1) 
fusion rules for the gaussian field with the fusion rules in the corresponding boundary CFT. 

2. In the ADE rational string theories, the matrix model realization of the boundary 
ground ring can be achieved by reformulating the corresponding ADE matrix chain in 
terms of A-type and S-type complex matrices. 
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Consider for example the simplest case of a theory of the A-series with central charge 
c = l~ p(p^i) ; which corresponds to 6^ = The corresponding string theory is described 
by the critical point of the Ap matrix model or by the tricritical point of the Ap-i matrix 
model. Both models can be incorporated in a single matrix model of the type (5.5) by 
replacing the Dynkin graph of A^o with that of Ap. The new model is defined in terms 
of p — 1 A-matrices labeled by the links of the Dynkin graph and p S-matrices labeled by 
the points of the graph. 



3. It is interesting that Dijkgraaf and Vafa []5^ proposed a similar matrix representation 
for the chiral ring for the c = 1 string theory compactified at the self-dual radius. Since the 



boundary operators can be considered as operators in a chiral CFT [|59| , |60[| , there might 



be a connection between the Witten's chiral ring and the boundary ground ring. 

4. The action of the boundary ground ring closes on the discrete subset s = + inb, 



m^n E This subset is relevant for the ZZ brane boundary states as discussed in 
||39| , p2| , |61| . [T0|] . One might ask oneself whether in the ZZ brane there is an analog of the 
boundary correlation functions in the FZZT brane. 
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Appendix A. Bulk and boundary reflection amplitudes 

The closed and open string tachyons of same momenta but opposite chiralities are related 
by the bulk and Liouville boundary reflection amplitudes [1^ 



Vp(+) = S^p+-%-\ Vp(-) = S^p-+^V^p+^ (A.l) 

[B^+Y^- =D^+-\si,S2) [B'pY''' =D'f+\si,S2) [B^iY^'. (A.2) 
In the normalizations (2.8) and (2.9) 

4+-^ = -1a^^ S^p-+^ = -b' a-""/' (A.3) 



D 



(+- 



(S1,S2) = - 



1 A^/' 



: S1+S2 



: Si— £2 



bS,{2P + b) s,(^q-P-^^'^ S,l^^-P 



A S1-S2 



(A.4) 
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where Ti,{x) is the Double Gamma function (see [|T^) and Sb{x) = Ti,{x)/Ti,{Q — x). The 
amphtude -D^ ^ = {Dp has the symmetries 



D^r-D^r- = d',-*'d%*' = t,'"^^ (A.5) 

^ sm2TTbP ^ sin^P 



and 

I^i7^(si,S2) = ^i.+"^(si,S2), (A.6) 

where the function \si,S2) is defined by (A. 4) with b replaced by 6 = 1/6. The 

reflection amplitudes for the degenerate momenta Pmn = \{rn/b — nb) are given by rational 
functions of z = Mcosh(7r6s) and z = M^/^ cosh(7rs/6). For example, 

Do-+ = 6^(6) = 62 D+- = lSil/b) = ^ 

Db/~2 = ih^(^l+^2) DZ+/2f, = sin(t%) (^l +^2) (A.7) 

7-) — I- _ sin(7rb^) gj-gg nH — _ 21-^2 

-1/26+6/2 sin(7r/fe^) 21-22 "^1/26-6/2 21-22" 



Appendix B. Recurrence relations for the closed string tachyon amplitudes 



Using the ring relations (3.2) and (3.6) one can obtain a set of recurrence equations for 
the correlation functions of the bulk tachyons 



n-1 



.k=l 



ccvi,:)(o)ccv),:;^(i 



(+) 



n+m — 1 

n 

j=n+2 



V 



(+) 



ccvL+^ (oo) 



(B.l) 



which generalize the recurrence relations for the resonance amplitudes obtained in |TI 
Here we assume that the neutrality condition is satisfied 



n+m 



(^0 - Pk 



2er 



(B.2) 



k=i 



The auxiliary function 

F(^Xj x\ P\ , . . . , Pn I Pn-\-l 1 •••I Pn-\-r 



n-1 



.fc=l 



n V 



Pk 



ccV^pJ+,iO) a_(x,x) ccV)r;^(l 



(+) 



m+n — 1 

n 

j=n+2 



V 



(+) 



ccV 



(+) 



,(oo) 



(B.3) 

does not depend on x and x. This can be proved by using d^a- = {Qbrst, b_ia_} and 
deforming the contour, commuting Qbrst with the other operators in (B.3) |Tl[]. Therefore 
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one can evaluate this function at x = and a: = 1 by using the fusion rules (3.2) or (3.5) 
and (3.6). As a result one obtains the recurrence relation 



G{Pi, Pn\Pn+l, Pn+m) — ^ G{Pi^ P„ + 6|P^_|_i — 6, P^+r 



m — 1 



— G{Pl, ...,Pn + b\Pn+2, Pn+j + Pn+1 — j;, 

Similarly, by inserting a+ we get the dual recurrence relation 

G{Pij P„|P^_|_i, Pn-\-m) = A G{Pi, Pn + j^\Pn+l — ^, Pn+m)- 
n-1 

— G{Pi, Pk + Pn + b, Pn-l\Pn+l, Pn+m) ■ 



k=l 



(B.4) 



(B.5) 



Note that the three-point function m + n = 3 coincides with the corresponding Liouville 
three-point function. The latter has been evaluated for generic momenta {i.e. without 
imposing the neutrality condition (B.2) by using a very similar argument in ||25|] . 



Appendix C. Recurrence relations for the boundary correlation functions with 
pure Neumann boundary conditions. 

The fusion rules above lead to a system of recursion relations for the disc correlation 
functions with arbitrary number of boundary states 



W{Pi, Pn\Pn+l, Pn+r 



'n-1 



.k=l 



n / 4:^ 



(+) 



n+m — l 



n B 



3=n+2 



'(+) 



c4+^ (oo) 



(C.l) 



Here three of the integrations are canceled by the volume of the global SL{2, IR) conformal 
symmetry of the upper half plane. Consider the auxiliary function 



-P ( X I Pi , . . . , Pn I Pn+1 5 • • ■ 5 Pn+r 



n-1 



.k=l 



cB^-\,{0)A.{x) cB'pZ'il 



(+) 



Pn+ 2 



n+m — l 



n « 



j=n+2 



?(+) 



cB^j^^ (oo) 



(C.2) 



Since 



dxF = {Qbrst, b_ia_} 
27 



(C.3) 



the function F does not actually depend on x and one can calculate it in two different ways 



by taking the limits x ^ and a; ^ 1 and using (4.3) and (4.6) (see [|TT| for a discussion 



concerning possible boundary terms). Note that the boundary operators are ordered and 
the fusion can be ] 
recurrence relation 



the fusion can be performed only with the operators Bp ■* and ^ . This leads to the 



W^D DID D ^ Vr(Pl,...,Pn+ +^'n+2 - ^,-Pn+3,---,-Pn+m) 

W (^-Tl, J-Vi|J-Vi_|_i, rn+m) — 



sin(27r6Pn+i) 

(C.4) 

or, after shifting the momenta, 

W/'D D fciD I 1 D \ _ ^(-Pl, -Pnl-Pn+1 + -Pn+2, -Pn+3, -Pn+m) 
VV , I^n — 21 "T 2& ' " " " ' -^ri+m ) — 



sin(27r6P„+i) 

(C.5) 

and similarly for which can be easily solved: 

n{n — 1) . m( m — 1) 

nr="i sm27r6(Pi + ... + P,)nr="i sinf (Ki + ... + Kz) 

Appendix D. Evaluation of the basic boundary correlation functions on the 
FZZT brane 

D. 1 . Three-point function 

For the three-point function 

the functional equations (4.21) and (4.22) read (with the assumption that in the first 
equation -Pi < — -P2 > ^ and in the seccond equation Pi < — |, P2 > |) 



sin (27r6P2) Wp^^p^,p^{s' , s ± ih, s") 



2M cosh 



s-\-s 



±iP') 



cosh 



7r6 



W 



Pl+2' 



p(s',s,s") (D.2) 



sin(^P2) Wp,,p,,p,{s'.s±ls") 
= 2M^/^' cosh 



±zP, 



cosh 



TT I S — S 



±iPn 



Wp , 1 1 pAs\ s, s") 



(D.3) 
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Their common solution is given by the boundary three-point function in pure LiouviUe 
theory (with the restriction P1-I-P2+-P3 = eo/2 on the momenta and rescaled to match with 
our definition of the boundary cosmological constant), which has been found exphcitly as 
an integral of products of double sine functions . If the three momenta are degenerate, 
that is of the form [r /b — sb) /2, the answer is a rational function of trigonometric functions 
12^ 



D.2. Two-point function 

The functional equations (4.21) and (4.22) are derived for the correlation functions of three 
or more operators. The two- and one-point functions in 2D QG are more subtle quantities 
whose correct normalization requires to tame the residual global conformal symmetry. To 
avoid facing this difficult technical problem, we will apply a trick: we will use the relation 
between the two-point function 



disc 



:D.4) 



and the three-point function with —Pi = P2 = P and P3 = |eo, and then apply (4.21) 
and (4.22). The relation is 



W 



Dp{s,s')-Dp{s,s") 



z{s') - z{s") 

In particular, taking s' = s" we have the more familiar relation 

dDp{s,s') 



dz{s') 



W_ppi^(s',s,s'). 



(D.5) 



(D.6) 



In this case eqn. (D.3) reads 

sin (^P) [Dp{s ± ib, s') - Dp{s ± ib, s")] 



2M cosh 



cosh 



[D t{s.s')-D^ tis^s' 
^ 2 



P-2 



where it is assumed that P > |, so that the shifted momentum has the same sign. This 
equation implies the following linear difference equation for two-point function: 



Dp{s±ib, s') 



2M cosh 



7r6 



s-\-s 



±iP 



cosh 



Tib 



±iP 



D 



p-b 
^ 2 



sin (27r6P) 



:D.7) 



The easiest way we know to prove that is to perform a Laplace transformation from the 
boundary cosmological constant z{s) to the physical length i = J dxe^^^^'^ for each segment of the 
boundary. 
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Similarly one obtains, for P > ^i the dual equation 



Dp{s±ls') 2MVf%osli 



cosh 



sm 



(D.8) 



The form of equations (D.7) and (D.8) suggests to look for a solution in a factorized form 



Dp{s, s') = bp{s + s')bp{s - s') 



(D.9) 



where Dp[s) should satisfy 



bp{s±ih) ^ 2^1/2 cosh[7rb(f ±zP)] 

•v/2sin(27r6P) 



^ 2 



(P>|). 



(D.IO) 



and 



D 



'2b 



A/2sin(2fP) 



(D.ll) 



In addition the two-point function should be real even function of P, s and s' . The solution 
is proportional to the Liouville reflection amplitude and is unique up to normalization: 



Dp{si,S2) 



1/b 



sin(27rP/6)|^l^l ^^^''^'^ ~ |sin(27rP6)| 



D^-t}{si,S2). (D.12) 



Here the normalization is fixed by the requirement that the amplitude is self-dual. 

D.3. One-point function 

One can proceed in the similar way to evaluate the boundary one-point function, or the 
disc loop amplitude 



using its relation with the boundary two-point function 

z{Si) - Z{S2) 

The explicit expression of -Deo/2 follows from the last formula in (A. 7): 

^1/6^-1 cosh(f s) - cosh(f s') 



(D.14) 



P>eo/2{s,s') = -b 



sin(|-eo) cosh(7r6s) — cosh(7r6s') 



(D.15) 
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Comparing (D.15) and (D.14) we get 



W{s) = b 



Mi/^'cosh(fs) 
sin(^eo) 



Similarly one obtains the dual loop amplitude 

^ 1 cosh(^_ 

sm(7roeo) 

D.4- The disc partition function of the FZZT brane 

The one-point function is related to the disc partition function S{z) by 

dS{A,s) 



W{s) 



dz 



z{s) = M coshTrbs. 



Integrating in z{s) gives for the disc partition function 
^2^^^1+1/6^ /cosh(7r(5s) cosh(7reos) 



^ ' 2sin(feo) 

sin(f eo) Qcq 
The dual partition function is 



Q eo 
(6cosh(7r6s) cos(|-s) — | sinh(7r6s) sinh(|-s)) . 



, fc-^M^+V''^ /cosh(7rQs) cosh(7reos) 



2sin(7r6eo) 



\ Q eo 

(— cosh(7r6s) cos(|-s) + 6sinh(7r6s) sinh(^s)) 



sin(7r6eo) Qbq 

The two partition functions are related by a Legendre transformation 
\ sin(|eo)-S'(s) + 6sin(7r6eo) S{s) + z{s)z{s) = 0. 



(D.16) 



(D.17) 



(D.18) 



(D.19) 



(D.20) 



(D.21) 



Appendix E. Bilinear equation for the loop amplitude 

The Virasoro conditions on the loop amplitude (5.16) are 



W^iz)+(b ^ 



z' W,{z)-W^{z') 



z — z' 



2T + W^^+i(2T-S') + W^x-i(2T-5') — = (E.l) 



dz' 



where the contour C\ encircles the support of the eigenvalue density [0, a] but not the 
interval [2T — a, 2T]. For amplitudes of the form (5.18) this equation becomes 



W\z)^ 6 - 



z' W{z)-W{z') 



2T + 2cos(7r62) W^(2r - 5') 
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dz' 
27H 



= 0. 



(E.2) 



The function W{z) is meromorphic with a cut along the support [0, a] of the eigenvalue 
distribution. Equation (E.2) implies the boundary condition along the cut 

W{z + iO) + W{z-iO)-2cos{nb^)W{2T-2)=2T (0 < z < d). (E.3) 

The integral equation (E.2) can be turned into a functional equation by the following 
trick ( sect. 3.2). Replacing z by 2T — z in (E.2) and adding the two equations one 
can turn the integral along the contour Ci into an integral along a contour Coo, which 
encircles both intervals [0, a] and [2T — a, 2T] and therefore can be expanded to infinity. 
Now the integration can be performed by the residue formula using the asymptotics 

W{z) = N/z + Wi/z^ + ... (z^oo) (E.4) 

which yields the the functional identity: 

W^{z) + W^{2T -z)-2 cos{Tih'^)W{z)W{2T - z) 



2T 



W{z) + W{2T - z) 



4T^C (E-5) 



z{2T-z) 

Here the constant C = C{N, T) is given by the integral 

T Jc^ 2m 

The critical value of C can be found by solving (E.5) for z = T, where it becomes 
algebraic. This gives 



^(^) = TT-^TT^ - T-?^' Cc = --^-T^. (E.7) 



2sin^(7r62) sin(7r&2) ' 4sin'(7r62) 
At the critical point the solution of (E.5) is given by 



~ rr. T , ^sinh(l -62)^^ 

z-T= — — , W = A + B ^— E.8 

cosh u smh u 

with A = o ■ 1,2 /o^ and B = . . .^/J ■ ^ .^n . The coefficient B is such that W(z) = 
ai z — oo. For given T the critical value of is determined by comparing the large z 
asymptotics of the solution (E.8) with (E.4): 

Nc = T'. (E.9) 

For A = N — Nc << Nc the constant C behaves as 

C = Cc + ciA + C2A'/'"t^-'/'" + ... 
and as ^2 < 1^ only the first term survives in the scaling limit A/T'^ — > 0. 
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